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Abstract 

The delta function gauge of the electromagnetic potential allows semiclassi- 
cal formulas to be obtained for the probability of exciting a single electron 
out of the ground state in an ultrarelativistic heavy ion reaction. Exact for- 
mulas have been obtained in the limits of zero impact parameter and large, 
perturbative, impact parameter. The perturbative impact parameter result 
can be exploited to obtain a semi-empirical cross section formula of the form 
a = A\n^+B for single electron ionization. A and B can be evaluated for any 
combination of target and projectile, and the resulting simple formula is good 
at all ultrarelativistic energies. The analytical form of A and B elucidates 
a result previously found in numerical calculations: scaled ionization cross 
sections decrease with increasing charge of the nucleus being ionized. The 
cross section values obtained from the present formula are in good agreement 
with recent CERN SPS data from a Pb beam on various nuclear targets. 
PACS: 34.90. -Kq, 25.75.-q 

I. INTRODUCTION 



In a recent work [|T| ionization cross sections were calculalated for a number of represen- 
tative cases of collisions involving ultrarelativistic Pb, Zr, Ca, Ne and H ions. The method 
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of calculation (on a computer) involved an exact semiclassical solution of the Dirac equation 
in the ultrarelativistic limit 0. A single electron was taken to be bound to one nucleus with 
the other nucleus completely stripped. The probability that the electron would be ionized 
in the collision was calculated as a function of impact parameter, and cross sections were 
then constructed by the usual integration of the probabilities over the impact parameter. 
The results of the probability calculations were used to construct cross sections for various 
ion-ion collision combinations in the form 



where A and B are constants for a given ion-ion pair and 7 (= 1/^/1 — v"^) is the relativistic 
factor one of the ions seen from the rest frame of the other. 

In Section II of this paper analytic results are derived for the probability that a single 
ground state electron will be excited in an ultrarelativistic heavy ion reaction. Exact semi- 
classical formulas are presented for the limits of zero impact parameter and perturbational 
impact parameters. In Section III the perturbational impact parameter analytical form is 
used as a basis to construct semi-empirical formulas for A and B. These formulas reproduce 
the previous numerical results for single particle ionization, and they illuminate the system- 
atic behavior of A and B with changing target and projectile ion species. Ionization cross 
sections calculated with Eq.(l) are then compared with data. 

II. IMPACT PARAMETER DEPENDENT PROBABILITIES 

If one works in the appropriate gauge 0, then the Coulomb potential produced by an 
ultrarelativistic particle (such as a heavy ion) in uniform motion can be expressed in the 
following form 0] 



b is the impact parameter, perpendicular to the 2;-axis along which the ion travels, p, z, and 
t are the coordinates of the potential relative to a fixed target (or ion), is the Dirac matrix, 
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a = Aln.'y + B 



(1) 





a the fine structure constant, and Z\ and v the charge and velocity of the moving ion. This 
is the physically relevant ultrarelativistic potential since it was obtained by ignoring terms 
in (b — p) [0 . Its multipole expansion is 

V{p,z,t) = aZr{\-OL:)b{z-t) 
^ ^ 2 cosm(/) 
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For b» p 



V{p,z,t) = 5{z - t)aZi{l - a,)2j- cos (p. (4) 

As will be shown in Section III, when b becomes large enough that expression Eq.(4) is inac- 
curate for use in calculating a probability, we match onto a Weizsacker- Williams expression 
which is valid for large b. Note that the b^ in the denominator of the logarithm in Eq. (2) is 
removable by a gauge transformation, and we retain the option of keeping or removing it as 
convenient. 

It was shown in Ref. that the 6 function allows the Dirac equation to be solved exactly 
at the point of interaction, z = t. Exact amplitudes then take the form 



a}(t = oo) = 6fj+ / dte'^^f-''^^'{(f)f\6{z - - a,) 

J—oo 

^^g-*aZiln{b-p)2 _ ^^1^.^ (5) 

where j is the initial state and / the final state. This amplitude is in the same form as the 
perturbation theory amplitude, but with an effective potential to represent all the higher 
order effects exactly, 

V{p, z,t) = -t6{z -t){l- a,)(e-^°^ii'^(''-^)' _ 



in place of the potential of Eq.(2). 

Since an exact solution must be unitary, the ionization probability (the sum of probabil- 
ities of excitation from the single bound electron to particular continuum states) is equal to 
the deficit of the final bound state electron population 

Em = i-Em (7) 

ion bound 

The sum of bound state probabilities includes the probability that the electron remains in 
the ground state plus the sum of probabilities that it ends up in an excited bound state. 
From Eq. (5) one may obtain in simple form the exact survival probability of an initial state 

^,W = 1(0.1(1 -«.)e-^"^^^°^'-^^'l0.)r- (8) 
By symmetry the az term falls out and we are left with 

P,(6) = |(0,|e-^"^^i'^(^-^)'|0,)r. (9) 
The ground state wave function (pj is the usual K shell Dirac spinor |^ 



with upper and lower components wave functions g and / 

72-1 -aZ-ir 



(10) 



g{r) = N^l + 72 r^^"^ e 



/(r) = -NJi - 72 r''^-^ e""^^" (11) 



where Z2, is the charge of the nucleus that the electron is bound to, 72 = y 1 — a^^'l, and 

2 (2aZ2)2^2+i , , 

= S — -11 -. 12) 

2r(72 + i) ^ ^ 

Let us first consider 6 = 0. We have 

p.(b = 0) = K^ile-^^"^^'"''!^,)^ = |(0,|e-2'"^^(''^"+''^('^'"^»|0,-)|2. (13) 

Putting in the explicit form of the upper and lower components for the K shell lowest bound 
state Dirac wave function and carrying out the integration we have 



or 



^.(& = 0) = ^ 



P,(6 = 0) 



r(272 + 1 - 2iaZi)T{l - iaZ^ ^ 



r(272 + i)r(|-mZi) 



(14) 



TTaZi ctnh(7raZi) 


r(272 + 1 - 2iaZi) 




r(272 + 1) 



(15) 



It is interesting to compare this result with a previous calculation of the probability of 
ionization in "close collisions" by Bertulani and Baur . For a one electron atom they find 

Pion{h < Xc/aZ2) = 1.8 a^Zl (16) 

where Xc = h/ rrieC is the electron Compton wavelength. If we take the low Zi limit of our 
expression Eq.(15) and then subtract it from one we obtain 

P.on(6 = 0) = (y - = 2-29 a^Z! (17) 

However our expression Eq.(15) only gives the flux lost from the initial state; some of 
that flux goes into excited bound states and is not ionized. From our previous numerical 
calculations we find that the actual ionization probabilities obtained either by summing up 
final continuum states or else by subtracting all the final bound states from unity were 76% 
- 80% respectively of the flux lost from the initial state. Thus if we multiply the constant 
in Eq.(17) by such a percentage we are in remarkable agreement with Bertulani and Baur 
for the perturbative limit. 

Now let us consider the case of 6 >> p. From Eq.(4) and Eq.(9) we have 

P,(6) = |(0,|e-2^"^^=°^('^)(''/^)|0,)|2. (18) 

Expanding the exponential up to we have 
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p.(b) = |(0j|l-2iaZicos(0)^-2a2Zi2cos2(0)^|0j)P (19) 
The term in cos{(f)) vanishes by symmetry, and integrating, we obtain 

5 



by ignoring the term in 1/6^. 

Both hmits, Eq.(15) for 6 = and Eq.(20) for b » p, are relativistically correct and 
thus correct for all Zi and Z2 since exact Dirac wave functions were used. 



III. A SEMI-EMPIRICAL FORMULA FOR SINGLE ELECTRON IONIZATION 

It is well known that the cross section for ionization of any pair of projectile and target 
species can be expressed as a sum of a constant term and a term going as the log of the 
relativistic 7 of the beam as seen in the target rest frame [0] 

aion = Alwy + B. (21) 

The cross section of this form is constructed from an impact parameter integral 

/ P(6),o„ h db (22) 

where P{b) is the probability of ionization at a given impact parameter. If all the flux lost 
from the initial state went into the continuum then Eq.(20) would provide the ionization 
probability at moderately large b 

We will take this form as a physical basis to build a semi-empirical formula for ionization. 
In any case we need to integrate the probability up to a natural energy cutoff. In order to 
do this we match the delta function solution Eq.(23) at some moderately large b onto the 
known Weizsacker- Williams probability for larger b by noting that if 6u; << 7 then 

and we can rewrite Eq.(24) in the Weizsacker- Williams form for large b 

Z?(l + 372 + 27|)A^.;2 
Pion{b) = 2-2 -^Ki (— ). (25) 

To perform the large b cutoff recall that to high degree of accuracy 
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/""^ ,.2/^&N, n 1 .O.68I7, , , ,0.681, 

— / Kl{—)h db = ln(— = 1117 + ln(— -). 26 

7^ Jbo 7 uJOo ujoq 

We immediately obtain the following expression for A 

^ = ^^(1 + 372 + 27!), (27) 

where A^, the square of the electron Compton wave length, is 1491 barns. However, as it turns 
out, uniformly for all species of heavy ion reactions, at perturbational impact parameters a 
little over 70% of the flux lost from the initial state goes into excited bound states and does 
not contribute to ionization. But since the ratio of flux going into continuum states to the 
total flux lost is so uniform we can use a flt to previously published numerical results ^ to 
obtain a semi-analytical form for A: 

A = (0.2869)^^(1 + 372 + 272'), (28) 

or in barns 
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A = 1792||(l + 372 + 272'). (29) 
^2 

Now one can use the second term in Eq.(26) to obtain a provisional expression for B 

B = Aln{——). (30) 

UJOq 

Obviously we need to evaluate u and to discuss 60 • ^ can be taken as the minimum ionization 
energy, 1 — 72, times a constant a little larger than one. One next observes that if Pionib) 
varies as 1/6^ the impact parameter integral has to be cut off on the low side at some value 60 
to avoid divergence. In fact the 1/6^ dependence continues down to the surface of the atom 
where other terms evident in Eq.(3) begin to contribute. The atomic size is just the electron 
Compton wave length divided by 0:^2 • In this region Pion{b) flrst rises faster than 1/6^ and 
then levels off to approach a constant change with 6 at 6 = [0]. One could try to add a 
low impact paramenter contribution to A based on Eq.(15) to our provisional form Eq.(30), 
but that turns out to unduly complicate things without improving the phenomenology. Our 
approach will be to set bo to an empirical constant divided by aZ2 



Eq.(30) now takes the form 



B = A\ii{ 



CaZ2 
1-72 



(31) 



Putting in two analytical fine tuning factors and fitting the remaining constant to the nu- 
merical results of Ref. we obtain a semi-analytical form for B: 



Table I expands a corresponding table from Ref. by adding cross sections of symmetric 
ion-ion pairs calculated with the formulas for A and B. There is good agreement between the 
formula values for the cross sections (first rows) and the numerical cross sections calculated 
by subtracting the bound state probabilities from unity (second rows) or calculated by 
summing continuum electron final states (third rows). For both A and B the agreement is 
also good with the Anholt and Becker calculations in the literature for the lighter ion 
species. However with increasing mass of the ions the perturbative energy dependent term 
A decreases in the formula calculations and in our previous numerical calculations, whereas 
it increases in the Anholt and Becker calculations. The greatest discrepancy is for Pb + Pb, 
with Anholt and Becker being about 60% higher. The reason that the A should decrease 
with increasing mass (actually Z) of the ions is explained by the 



factor in the formula for A (and thereby B also). As we noted before, perhaps the discrepancy 
between our A decreasing with Z and the Anholt and Becker A increasing with Z is due 
to the fact that Anholt and Becker use approximate relativistic bound state wave functions 
and the present calculations utilize exact Dirac wave functions for the bound states. For 
the term B (which has the non-perturbative component) the agreement is relatively good 
between all the calculations. 

In the perturbative limit (small Zi, Z2) the cross section formula goes over to 



B = A-i\'^\l - a'Z^Z^f'Hxii 



2.37aZ2 

1-72 



(32) 




(33) 



a = (0.2869)87rA^-|ln 



Z^ 



2.377 

aZ2 




In 



2.377 

aZ2 



(34) 
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By way of comparison, Bertulani and Baur |^ using the equivalent photon method and 
taking the contribution of 6 > Ac/ found 



a = 4.9Ac^||ln-^ (35) 

Z9 a/j' 



for this case of ionization of a single electron. 

Table II shows results of the calculation of B (multiplied by Z\jZ'^ for a number of 
representative non-symmetric ion- ion pairs. (Since A is perturbative, scaling as its 
value can be taken from Table I for the various pairs here.) Once again there is good 
agreement between the formula values for the cross sections (first rows) and the numerical 
cross sections calculated by subtracting the bound state probabilities from unity (second 
rows) or calculated by summing continuum electron final states (third rows). The only 
notable disagreement is with Anholt and Becker for Pb targets. 

The availabilty of the present semi-empirical formula facilitates a comparison with avail- 
able CERN SPS data. Calculations with the formula are in considerably better agreement 
with the data of Krause et al. |^ for a Pb beam on various targets than are the Anholt and 
Becker numbers with or without screening. Note that in this case the role of target and 
beam are reversed. It is the single electron Pb ion in the beam that is ionized by the various 
nuclei in the fixed targets. The formula numbers do not include screening, which should in 
principle be included for a fixed target case. However, one might infer from the Anholt and 
Becker calculations that the effect of screening is smaller than the error induced by using 
an approximate rather than proper relativistic wave function for the electron bound in Pb. 

Note that the formula has not been fit to experimental data. It is compared with ex- 
perimental data. The "empirical" aspect of this formula refers to adjusting the formula to 
previous numerical calculations of Ref. |1| 

At RHIC the relativistic 7 of one ion seen in the rest frame of the other is 23,000, and 
of course there is no screening, so the present formula should be completely applicable. The 
present formula predicts a single electron ionization cross section of 101 kilobarns for Au + 
Au at RHIC. The corresponding cross section from Anholt and Becker is 150 kilobarns. 
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TABLES 



TABLE I. Calculated Ionization Cross Sections Expressed in the Form A\n^ + B (in barns) 





Pb + Pb 


Zr + Zr 


Ca + Ca 


Ne + Ne 


H + H 


A 


Formula 


8400 


10,212 


10,618 


10,718 


10,752 






8680 


10,240 


10,620 


10,730 


10,770 






8450 


9970 


10,340 


10,440 


10,480 




Anholt & Becker 


13,800 


11,600 


10,800 


10,600 


10,540 


B 


Formula 


14,133 


27,375 


36,623 


44,638 


69,629 




1 — Hbnd e 


14,190 


28,450 


38,010 


46,080 


71,090 






12,920 


27,110 


36,530 


44,430 


68,780 




Anholt & Becker 


13,000 


27,800 


37,400 


45,400 


70,000 
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TABLE II. Calculated values of the scaled quantity {Z^jZ'^B for non-symmetric combinations 
of colliding particles. The second nucleus (Z2) is taken to be the one with the single electron to 
be ionized. Since Anholt and Becker cross sections without screening are completely perturba- 
tive, their values of of B also can be taken from Table IV, and are repeated here for convenient 
comparison. 





H + Ne 


H + Ca 


Ca + H 


H + Zr 


H + Pb 


Pb + H 


i VJllilLllCl 


tit: ^ \ XU 


36 8Q0 


69 462 


28 226 


1 6 487 


66 539 


1 ~ Tlihnd ^ 


46,150 


38,270 


70,820 


29,440 


17,090 


67,550 




44,490 


36,790 


68,520 


28,070 


15,680 


65,330 


Anholt & Becker 


45,400 


37,400 


70,000 


27,800 


13,000 


70,000 




Pb +Ne 


Ne + Pb 


Pb + Ca 


Ca + Pb 


Pb + Zr 


Zr + Pb 


Formula 


42,308 


16,313 


34,503 


16,097 


25,751 


15,592 


1 ~ Hhnd 6 


42,560 


17,030 


34,720 


16,870 


26,010 


16,250 




41,000 


15,690 


33,330 


15,530 


24,730 


14,930 


Anholt &; Becker 


45,400 


13,000 


37,400 


13,000 


27,800 


13,000 
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TABLE III. Cross sections for the ionization of a 160 GeV/A one electron Pb projectile {Z2) 
by various fixed nuclear targets (Zi). Unlike in Table II, here the appropriate factor has 

been included. Cross sections are given in kilobarns to match the format of the CERN SPS data. 



Target 


^1 


Formula 


Dr^D JJata 


Anholt & Becker 
(with screening) 


Anholt & Becker 
(no screening) 


Be 


4 


0.14 


0.14 


0.24 


0.20 


C 


6 


0.32 


0.31 


0.49 


0.45 


Al 


13 


1.5 


1.3 


2.0 


2.1 


Cu 


29 


7.4 


6.9 


9.0 


10.5 


Sn 


50 


22 


15 


25 


31 


Au 


79 


53 


42 


60 


78 
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